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Low PMEPR OFDM radar waveform design using 
the iterative least squares algorithm 

Tianyao Huang and Tong Zhao 


Abstract —This letter considers waveform design of orthogonal 
frequency division multiplexing (OFDM) signal for radar appli¬ 
cations, and aims at mitigating the envelope fluctuation in OFDM. 
A novel method is proposed to reduce the peak-to-mean envelope 
power ratio (PMEPR), which is commonly used to evaluate the 
fluctuation. The proposed method is based on the tone reservation 
approach, in which some hits or suhcarriers of OFDM are 
allocated for decreasing PMEPR. We introduce the coefficient of 
variation of envelopes (CVE) as the cost function for waveform 
optimization, and develop an iterative least squares algorithm. 
Minimizing CVE leads to distinct PMEPR reduction, and it is 
guaranteed that the cost function monotonically decreases hy 
applying the iterative algorithm. Simulations demonstrate that 
the envelope is significantly smoothed hy the proposed method. 

1. Introduction 

As a multi-carrier modulation waveform, OFDM is em¬ 
ployed both in communication and radar applications. It has 
received growing attention from radar community for its 
excellent performance since the pioneer work on multi-carrier 
phase-coded (MCPC) signal by Levanon [2] (Chapter 11). 
OFDM is regarded as a low interception probability (LPI) 
waveform [3] (Chapter 5) for its instantaneous wideband prop¬ 
erty. The frequency diversity in OFDM provides additional 
information on observed targets, and thus increases detec¬ 
tion probabilities of the targets [4]. The wideband waveform 
mitigates possible fading and resolves the different multipath 
reflections [5, 6]. In addition, due to advantages both in radar 
and communication applications, OFDM is also investigated 
for joint radar-communication systems [7, 8]. 

A severe drawback of OFDM is the envelope fluctua¬ 
tion problem. Superposition of different subcarriers generally 
brings notable variations in amplitude. The variation in instan¬ 
taneous power leads to a large power backoff and reduces the 
efficiency of the amplifier, or results in nonlinear distortion. 
This letter focuses on mitigating the envelope fluctuations. 

The peak-to-mean envelope power ratio (PMEPR) or peak- 
to-average power ratio (PAPR) is usually employed to measure 
envelope fluctuations. In order to alleviate power back off or 
signal deformation, many approaches have been proposed to 
reduce PMEPR/PAPR. 

Some certain amplitude-phase modulations to generate low- 
PMEPR/PAPR OEDM waveforms are investigated [2, 9, 10]. 
Golay complementary sequences yield a signal of PAPR 
bounded with 2 [9]. Newman and Narahashi-Nojima phasing 
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weights produce PMEPRs around 1.8 [2]. Schroeder method 
considers both amplitude and phase weight [2]. These so¬ 
lutions have closed-form expressions and are easy to apply, 
however, the resulting PMEPRs are not close enough to the 
achievable optimum [2]. In addition, since the modulations are 
immutable, the flexibility of these waveforms are restricted. 

Communication community has developed different ap¬ 
proaches for PAPR reduction, which allow arbitrary digital 
data to be transmitted. Comprehensive surveys can be found 
in [9, 11-13]. These algorithms are mainly divided into two 
categories: distortionless and distorted methods. Compared 
with distorted approaches like clipping and companding trans¬ 
form techniques, distortionless methods avoid additional in- 
band and out-of-band distortion to the signals, preserve the 
orthogonality between subcarriers, and make it possible to 
independently process each subcarrier of radar returns. 

Among these distortionless approaches, tone reservation 
(TR) algorithm draws special interests. In a TR algorithm, 
some bits or subchannels are used to carry arbitrary informa¬ 
tion, which provides flexibility for radar waveform selection, 
and the rest are reserved to smooth the fluctuated envelope. 

This letter presents a novel TR algorithm to reduce the 
PMEPR of an OEDM signal. The cost function and the 
solution are different from the standard TR method [9]. In 
the standard TR approach, only the peak of the envelope is 
minimized. Note that the envelope fluctuation relies on overall 
extension of the waveform rather than a single sample. We 
choose the variance of the envelope as the penalty function, 
which considers both peaks and valleys of the envelope. And 
an efficient iterative least squares algorithm is proposed cor¬ 
respondingly. Theoretical analyses and numerical experiments 
demonstrate the effectiveness of the proposed algorithm. 

The rest of this letter is structured as follows. In Section 
II, the signal model of OEDM waveform and the definitions 
of PAPR and PMEPR are introduced. Section III models 
PMEPR reduction as an optimization problem and an iterative 
algorithm is developed. In Section IV, simulation results 
are presented to examine the performance of the proposed 
algorithm. Section V briefly concludes this paper. 

II. OEDM Signal Model 

Consider a complex-valued OEDM waveform containing 
N carriers transmitted simultaneously with a frequency in¬ 
terval of A/. The frequency of the nth carrier /„ is nAf, 
n = 0,1,..., N —1. Each carrier contains M amplitude-phase 
modulated bits with a duration of tb- The bit duration 4 is 
set as 1/A/, which yields the orthogonality between bits in 
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different carriers. The pulse width Tp is Mtb- Denote the mth 
code in the nth carrier as an,m, m = 0,1,..., M — 1, and the 
superposition of all carriers is 

N-l M-1 

g{t) = X! X! On.mrect (f - mfb) (1) 

n—0 m—0 

where rect(-) is a rectangle function defined as 


is a Fourier transform matrix with Fkn = as its fcth 

row, nth column element. Let x = [atj, ..., € 

£OsNM q ^ ^ G then 

X = Aa, (7) 

where A = diag(F, F,...,F) G cO^nmxnm^ 
this paper is to design codes a which smooth \x\. 


rect(t) 


1,0 < f < 4, 

0, else. 


( 2 ) 


OFDM waveform suffers from a shortcoming that the am¬ 
plitude | 5 (f)| varies along with time. Note that in an actual 
radar system, real-valued signal g{t) = Re [(/(f)] is transmitted 
and amplitudes |p(f)| is of pratical interest. To enhance the 
detection range, power amplifier is always contained in the 
radar signal generator. High peaks encountered in |p(f)| could 
drive the power amplifier into saturation, which produces 
nonlinear distortion and interferences among subcarriers [9]. 
To avoid operating the power amplifier into saturation, input 
backoff (IBO) is applied [9]. However, IBO shortens the 
detection range since the average power of the signal is 
decreased. It is preferable to reduce the envelope variations 
of the waveform, which remains the signal power unchanged 
and relieves the signal distortion simultaneously. 

For a real-valued waveform p(f), the peak-to-average power 
ratio (PAPR), defined in (3), is usually adopted to characterize 
the amplitude fluctuations. 


PAPR^(() — 


maxt |g(f)p 


(3) 


The discrete-time PAPR is usually easier to handle than 
the continuous-time PAPR. Assume the sampling rate is 
fs = OsNAf. When Og > 4, the discrete-time PAPR well 
approximates the continuous-time counterpart [11]. 

For the convenience of using complex-valued waveform, 
PMEPR is selected instead of PAPR. And in low-bandwidth 
cases (i.e. the bandwidth is far less than the radio frequency), 
minimizing PMEPR is a good approximation to finding lower 
PAPR [2] (Chapter 11). Denote the sampled waveform as 
xi = g{l/fs), ( = 0,1,..., OsNM — 1, and the discrete-time 
PMEPR is expressed as 


PMEPR,^, 


max; |x;p 


(4) 


where Ei[\xi\‘^] = represents the 

mean of jx/p. Let I = k + mOgN, k = 0,1,..., OgN — 1, 
and xi can be rewritten as 


N-l 

Xl = Xk+mOsN = ^ (5) 

n—0 

Denote X-^i = [XmOsNiXi-\-mOsNj • ■ ■ iXOsN—l-\-mOsN] ^ 

'j’jjg signal within the mth bit can be presented by 


(6) 

where dm = [ao.m, oi.m, • • •, aAr-i,m]^ € corresponds 
to the modulation symbols in the mth bit, and F G (^^OsNxn 


HI. Problem Eormulation 


A TR algorithm divides a into two sub-vectors, dgi and 
dgR, containing informative symbols and reserved symbols, 
respectively. Sets and are disjoint subsets of = 
S^[JS^ = {0,1,NM — 1}. Symbols in dgi can be 
arbitrary. Given dgi, symbols in dgn are reserved to mitigate 
envelope fluctuations. 

Informative bits can be carefully designed for some special 
purpose. Eor example, in a joint radar-communication system 
these bits bear communication information. In other cases 
without specific constraints, dgi can also be randomly chosen 
to increase the diversity of the waveform. 

Note that in a radar, energy of both dgi and dgn are useful 
and important. The performance of target detection mainly 
depends on the total energy of the transmitted dgi and dgn. 
Echoes of both symbols are received and processed by the 
radar, and equally contribute to target detection. TR method 
allocates energy over carriers, but does not directly affect the 
total energy of transmission. 

A desirable criterion is to minimize PMEPR, i.e. 


max; \xi\^ 

■ EiWxiV] 


s.t. X — Agidgi -\~ -AgRdgi 


(8) 


where Ag denotes columns of A indexed in the set S. To 
abbreviate expressions, let j3 = Ei [|a;;|], c = Agidgi, b = 
dgR and B — AgR. Due to the denominator of the objective 
function, (8) is rather difficult to solve. 

A standard tone reservation approach [9] simplifies (8) by 
ignoring the denominator, which yields 


min||x||oo, s.t.x = c + Bb, (9) 

b 


where |jx||oo = max; |x;| is the infinity norm of x. Linear pro¬ 
gramming solves (9). A variation of (9) minimizes Ei [jx;]^] 
instead [14]. Both criteria seek to decrease the peak values in 
jxj. However, the valley values of jxj, which also contribute 
to envelop variations, are not increased by these two criteria. 
It is not guaranteed that a sequence b leading to a constant 
envelope is the optimizer of the criteria. 

This letter proposes a novel criterion which minimizes the 
variance of envelopes, and both peak and valley values are 
taken into account. We introduce coefficient of variation of 
envelope (CVE) to evaluate the fluctuation degree of x, which 
is defined as 


CVE,, = 


El 


(|a:;| -Ep[\xp\ 


{EiWxiWY 

CVE is closely related to PMEPR. It is satisfied that 


( 10 ) 


1 < ^PMEPR,, < + 1, (11) 
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where the equalities holds if and only if x has a constant 
envelope; see proofs in Appendix A. When the envelope is not 
constant and the number of samples O^NM are large, bounds 
in (11) could be loose. However, it still can be indicated from 
(11) that 1) codes in b which result in a constant envelope, if 
the codes exist, lead to the minimum of CVE, and 2) a small 
PMEPR is guaranteed if CVE is small. Omit the denominator 
in (10), optimization problem becomes 


min El 
b 



iEp[K\]f 


s.t. X = c + Bb, 


( 12 ) 


or equivalently, 

O^NM-l 

min E {\xi\- Ep[\xp\])'^ , s.t. X ^ c + Bb. (13) 
1=0 


It still holds that codes in b making the envelop constant are 
the optimal solution of (13). Reformulate (13) in a matrix form 

min {Bb + c- {Bb + c- /3e^®) , 

" 1 (14) 

s-t. P = = ABb + c), 

where || • || i denotes the norm of a vector, and Zx is the angle 
of complex-valued x. A closed-form solution to (14) is hard 
to obtain due to involved complex and nonlinear operations. 

Symbols in b and its functions /3, 9 are calculated alter¬ 
nately. Assume that at the ith iteration, is obtained, then 

+‘'Hi’= ABb^^^+c). (15) 

Substitute j3 = and 6 = in (15) into (14), and (14) 
reduces to a least squares problem, of which the solution is 

b(i+i) = __Bt , (16) 

where B^ is the pseudo inverse B^ = [B^B'^ ' B^. 

It is proved that the cost function in (14) is monotonically 
nonincreasing versus the iteration counter z; see Appendix B. 

Remarks. Partitioning of the set affects the achiev¬ 
able PMEPR level. Intuitively, increasing the number of free 
variables, i.e. \S^\, benefit lowering the PMEPR. But in 
some applications, choice of are rigid. Eor example, in 
a jamming environment, jammed subcarriers are contained in 
and the corresponding bits are set zeroes to evade jammer. 
Without such purposes, division of is arbitrary and has 
the potential to be optimized during the waveform design 
processing. We may remain further discussions as the future 
work. 


IV. Numerical Results 

Numerical experiments are performed to demonstrate the 
advantages of the proposed method over the standard TR 
approach (9) and its variation [14]. To solve problem (9) and 
its variation [14] we apply CVX, a package for specifying 
and solving convex programs [15, 16]. These three algorithms 
are denoted as ’TR-CVE’, ’TR-max’ and ’TR-Elxl"^’, corre¬ 
sponding to their cost functions, respectively. In the following 
experiments, the baseband bandwidth is B = 50 MHz, and 
the sampling rate fs = 500 MHz. 


A. Convergence 

Suppose that there are = 6 subcarriers, of which each 
contains M = 1 bit. Two subcarriers, n = 2,3, are predefined 
with 02,0 = 03,0 = 1- The rest subcarriers are reserved as free 
variables for PMEPR reduction. Eor the proposed method, set 
the initial value as b^^'> = 0, and alternately calculate /3, 6 
with (15) and the symbols b with (16), respectively. After 800 
iterations, the final symbols b are obtained. Envelopes of the 
optimized waveforms and the initial counterpart are depicted 
in Eig. 1. Note that envelopes are not normalized, and the 
average powers are not 1, ^ 1. 

As shown in Eig. 1. The fluctuations in the envelope declines 
significantly using the proposed method. PMEPR decreases 
from 2 to 1.05. The standard tone reservation approach and 
its variation decrease peak values in \xi\, but PMEPRs are 
slightly reduced to 1.30 and 1.51, respectively. 



Fig. 1. Envelopes of the initial and optimized waveforms. 


The roots of PMEPRs, the corresponding upper bound in 
(11) and the cost function of (14) versus the iteration number 
are calculated and shown in Pig. 2, which demonstrates the 
monotonicity of the objective function in (14). PMEPRs are 
not monotonically reduced, since PMEPR is slightly different 
to the cost function. And the upper bound is not tight (values 
above 1.5 are not displayed). However, PMEPR, its upper 
bound and the cost function have similar trend. When the 
bound and the cost function are small, PMEPRs are also small. 


B. CCDF of PMEPR 

In this subsection, the informative subcarriers and symbols 
therein are randomly selected, and the complementary cumu¬ 
lative distribution function (CCDP) of PMEPR is examined. 
We consider two cases N = 6 and N = 10, among which two 
are randomly selected as informative subcarriers. Set M = 10, 
and the codes an,m C 0 . 3 ’ ^re equally likely valued among 
{— 1, l,j, —j}. The CCDPs of PMEPRs are displayed in Pig. 
3, which is obtained with 5 x 10^ Monte-Carlo trials. Pig. 
3 shows that the proposed iterative algorithm leads to lower 
PMEPR levels than the existing approaches (9) and [14]. 
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Fig. 2. Roots of PMEPRs, the upper bounds in (11) and the cost function 
in (14) versus iteration number. 



Fig. 3. CCDFs of PMEPRs for random informative symbols agi. 


C. Ambiguity function and detection performance 

In this subsection, the ambiguity functions and target detec¬ 
tion performances of the optimized waveforms are discussed. 
There are = 6 subcarriers with M = 10 bits in each carrier. 
Two carriers, n = 2,3, are predefined, and Chu codes are 
chosen with its mth phase as an,m = , n = 2 ,3, where 

7=1 for n = 2 and 7 = — 1 for n = 3. The rest bits are 
reserved as free variables for TR algorithms. We also examine 
’Uniform OFDM’ with the rest bits of a uniform amplitude 
but random phases an,m = n = 0,1,4,5, where 

<pn,m is independently uniformly distributed over [0,27r). The 
ambiguity functions are shown in Fig. 4 ^ Fig. 7. 

As shown in these figures, waveform designed with TR- 
CVE algorithm has similar ambiguity function with that of 
’Uniform OFDM’ waveform. And waveforms by TR-max and 
TR-E|a;|^ have high sidelobes. 

We also examine the detection performance of these wave¬ 
forms. Waveforms are normalized such that the average power 
is unique. An ideal point target is considered. The noise power 
is set as tr^ = 1. Noisy radar echoes are processed with 
matched filter. The results are compared with thresholds. Eor 
different OEDM waveforms, thresholds are selected such that 
the false alarm is 10“®. The radar cross section of the target 



-10 -5 0 5 10 

Delay (4) 


Fig. 4. Ambiguity function of the OFDM waveform produced with TR-CVE 
algorithm. 



-10 -5 0 5 10 


Delay (tb) 

Fig. 5. Ambiguity function of the OFDM waveform produced with TR-max 
algorithm. 
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Delay (fj,) 

Fig. 6. Ambiguity function of the OFDM waveform produced with TR-E|fc|^ 
algorithm. 


varies to change the signal noise ratio (SNR). The detection 
probabilities of the target is shown in Fig. 8 . 

As shown in Fig. 8 , the detection probabilities of different 
waveforms are close to each other. 
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Delay (4) 

Fig. 7. Ambiguity function of the ’Uniform OFDM’ waveform. 



Fig. 8. Ambiguity function of the ’Uniform OFDM’ waveform. 


V. Conclusion 

This letter considers a tone-reservation based algorithm to 
lower the PMEPR of an OFDM waveform for radars. We 
propose coefficient of variation of envelopes (CVE) as a metric 
to measure the envelope fluctuation, and an iterative least 
squares method is devised consequently. Simulations show 
that the algorithm significantly improve the performance on 
PMEPR reduction. 


/VFULNUIX A. 


Since max/ \xi\^ > Ei \xif > {Ei [|x/|])^, it holds 


1 < \/PMEPR^, < max\xi\/Ep [|a:p|], (17) 

max; \xi\ _ ^ _ max; \xi \ - Ep [|xp|] 

Ep [I^CpI] Ep [|a;p|] 

< 

= CVE,,, 

and then we have (11). 


OsNM-l 

E 




(18) 


Appendix B 

Define the cost function in (14) as f{b,P,0). We prove 

/ , (19) 

where fel*! is an arbitrary initial value and the rest variables 
are obtained with (15) and (16). 

Proof: Since is a solution of least squares, it obeys, 

/ . ( 20 ) 

Then we need prove 

Define a cost function 


( 21 ) 


L i/3,6) = 1/2 {y - /3e^«) ^ (y - /3e^«) , (22) 

where y G is a nonzero vector, ,5 C M, 0 S R^, and each 
entry of 8 belongs to [0, 27r). Set y — + c. To prove 

(21), we only need to prove that the minimizer of (22) is (15). 

It can be found lim L {/3,0) — +oo, which indicates that 

^—¥00 

the global minimizer of (22) is an extreme point. Calculate 
the partial derivatives 


dL 


= NI3-Rt , 


^ = -/3Im (e © y) 


(23) 


(24) 


where © is element-wise multiplication. The Hessian matrix 


H = 


N —Im © y^j 

—Im (e”-’® © y) /)diag (Re © y)) 


■ (25) 


Solve = 0 and ^ = 0, and the Hessian matrix becomes 


de 


i?l||=o.|L=o = diag([W,±/3|yo|,---,±/3|yAr_i|]y (26) 
The Hessian matrix (26) is positive definite, if and only if 


P = ^\\y\\i,d = 


which yields H = diag(^ [V,/3|y|^] 
or (15), is the minimizer of (22). 


(27) 

and indicates that (27), 
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